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Abstract 

We analyze a simple textbook approach to nonlinear oscillators 
proposed recently, disclose its errors, limitations and misconceptions 
and complete the calculations that the authors failed to perform. 

1 Introduction 

In a recent article Ren and He [1] proposed a simple method for the ap- 
proximate calculation of the period of nonlinear oscillators and applied it 
to three rather trivial toy models. In this paper we discuss this remarkable 
scientific contribution that is another illustrative example of the new trend 
in mathematical physics promoted by certain journals. 
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2 Their method and our improvements 

Ren and He [1] chose the dimensionless equation of motion 

u" + fiu) = (1) 

with the initial conditions 

m(0) = A, u\0) = (2) 

where the prime indicates differentiation with respect to the independent 
variable t. 

The authors gave the following recipe: "We always choose cosine or sine 
function as a trial function for nonlinear oscillators. For the above initial 
conditions, we choose 

u = A cos ujt (3) 

where u is the angular frequency of the nonlinear oscillator to be further 
determined. Substituting Eq. ([3]) into Eq. ([1]) results in" 

u" + fiAcosut) = (4) 

It is unbelievable that the authors did not realize that this equation does 
not apply to an arbitrary nonlinear oscillator for all values of t. In fact it is 
suitable only for the Hooke's force f{u) = u'^u. The authors proceeded as 
follows: "Integrating Eq. (jl]) twice with respect to t, we have 

u' = - [ f{u)dt (5) 



and" 







u{t) = [ u'dt (6) 
Jo 



Another incredible mistake, this equation does not satisfy the first initial 
condition in Eq ([2]). 

In the discussion and conclusions Ren and He [1] admirably argued that 
"For an oscillator with initial conditions ([2]) we have" 



u{0) = A, u{T/A) = 0, u{T/2) = -A, u(3T/4) = 0, u{T) = A (7) 



However, they did not bother to make it explicit that this set of equations 
is valid only for a nonlinear oscillator with odd force f{—u) = —f{u). It is 
also probable that they were not aware of this obvious fact. Fortunately the 
authors only treated such particular cases. 

Before proceeding with the discussion of this admirable piece of work, we 
recall a well known result in classical dynamics. If we multiply Eq. ([T]) by u' 
and integrate we obtain the textbook expression 



where E is a constant of the motion and dV/du = f. 

To continue with the discussion of the paper by Ren and He [1] we derive 
the correct expression that in our opinion is the basis of their approach. If 
we integrate Eq. ([T]) twice and take into account the boundary conditions 
already indicated above we obtain 



that those authors supposedly tried to derive. The main idea behind their 
approach is that if one introduces an approximate solution into the right- 
hand side of Eq. (Q the result is expected to be an improvement. To obtain 
an approximate analytical expression for the period T the authors resorted to 



u 
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- + V{u) = E 



(8) 
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the condition u{T /4) = and in the conclusions they stated that "The sug- 
gested solution procedure is valid for conservation systems with unchanged 
amplitude". In fact, any equation of the form ([1]) represents a conservative 
system and exhibits a constant amplitude given by V{u) = E that deter- 
mines the turning points. Did the authors know it?. As indicated above 
u{T/A) = is valid only for symmetric problems V{—il) = V{u) or odd 
functions f{—u) = —f(u). In fact, they only considered examples with that 
property without stating it explicitly in the general presentation of the ap- 
proach. To facilitate the discussion below it is worth noticing that for an odd 
force A > u{t) > if < t < T/4. 



3 Their examples and our improvements 

The first example chosen by Ren and He [1] is the well known and widely 
studied Duffing oscillator: 

f(u) = u + eu^ (10) 

Undergraduate students learn how to obtain satisfactory approximate solu- 
tions to this equation in most textbooks on classical mechanics. As shown 
below, the results of Ren and He [1] are of such kind, except that they were 
published in a research journal and were derived in a sloppy way. Appar- 
ently, those authors did not realize that the model parameters e and A do not 
appear separately in the results but in the form of the only relevant quantity 
p = eA"^, an observation that greatly facilitates the discussion of the results. 
They first substituted Eq. with u = 271 /T into the right-hand side of 
Eq. ([9]) and obtained an improved trajectory as well as the approximate 
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period [1] 

rra(p) = -JL= (11) 

This expression provides a reasonable approach to the actual period for all 
values of p, and in particular an acceptable estimate of the limit 

lim VpT(p) = Too = 7.416298709 (12) 

p— >oo 

In fact, Eq. (HH) gives us t]^ ^ 7.12. 

How did they derive a reasonable result from the wrong equation ([6])? 
Simply by the addition of another error that corrects the first one. Notice 
that the left- and right-hand sides of equations (10) and (11) in their paper 
do not match. They added convenient integration constants in the last steps 
to satisfy the boundary conditions. 

In order to improve this first estimate Ren and He [1] then tried the 
ansatz 

u{t) = Ai cos{ujt) + A2 cos{3(jjt) (13) 

where Ai + A2 = A. The additional coefficient requires an additional condi- 
tion and the authors chose 

n"(0) + /(A)=0 (14) 

In this way they derived a complicated system of equations that suspiciously 
they did not try any further. Besides, at first glance their results do not 
appear to be functions of p alone. For this reason in what follows we derive 
a suitable expression and verify if it is more accurate than the first approxi- 
mation discussed above. Starting from the same premises our expression for 
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the period results to be 

125TV (p + 1)^ - reseTT^T^ (p + l)^ + UOtt^T^ (p + l) (1607p + 735) 
-647r^T2 (48851p + 55125) + 127008007r^ = (15) 

This equation has four real roots -T^^^^^{p) < -T^^^^^{p) < < T[2'i](p) < 
and T^'^'^^(p) is the desired approximation that satisfies T[^'^](0) = 2%. 
A straightforward calculation shows that the approximation to Too is a 
root of 

125T^ - 76567r2T^ + 192 8 407r^T^ - 31264647r''T^ + 127008007r^ = (16) 

We thus obtain t]^'^' ^ 7.44 that is in fact more accurate than T^. 

If the Duffing equation is a textbook exercise, the second example studied 
by Ren and He [1], given by 

f{u) = e sgn{u) (17) 

is ridiculously trivial because it is almost impossible to miss the exact result 
(remember that this paper is published in a research journal). We know that 
sgn{u) = 1 for < t < T/4; therefore, if we substitute f{u) = e into the 
right hand-side of Eq. we obtain u{t) = A — and the exact period 
T = 4^/2X76. Any trial function that is positive definite in this interval 
leads to the same result, even if it does not satisfy the initial condition. 
The third example is 

f (u) = Lulu + eu\u\ (18) 
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but it is sufficient to consider f{u) = UqU + ev? by virtue of the argument 
given above. Notice that in this case the independent model parameters are 
tuo and p = eA. If we choose the trial function ([3]) we obtain their result 



from which it follows that ^ 6.75. 

If we use the trial function (fT3l) we easily obtain and improved expression 
for the period 



{971^ - 16) +8n^T^ (l6 - 457r2) + 167r^T^ (36971^ + 1136) -2949127r^ = 

(21) 



In this case there are only two real roots — Tio and T<L ^ 6.867 and the 
agreement of the positive one with the exact result Too = 6.868663935 is 
remarkable. 

4 Conclusions 

The method proposed by Ren and He [1] could be thought of as an undergrad- 
uate exercise on classical mechanics with the limitation that it is valid only 
for odd forces because of the necessary condition u{T/4) = 0. Their paper 




(19) 



TV {oOq + 2ujIp + p2) (9vr2 - 16) 

-S'k'^T^ [2hQuj^ + Ihulp (371^ + 16) + p^ [Ahix'^ - 16)] 

+167r^T2 [5120cu^ + p (3697r2 + 1136)] - 2949127r^ = (20) 



and the corresponding limit p — > oo 
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may be of some pedagogical value as an exercise for undergraduate students 
who may be asked to find as many mistakes as possible. The lecturer may 
even organize a kind of competition among groups of students dedicated to 
such a task. We think that it can be a hilarious class. The second example in 
that paper is trivial and, therefore, only useful as a curiosity for beginners. 

We have carried out the second-order approximation for the Duffing os- 
cillator that those authors left unfinished, and also showed how to do that 
calculation for the third example that they never tried. There is no need to 
say that the results of that paper have no serious utility whatsoever for actual 
research in the field of nonlinear oscillations. However, it is not surprising 
that such a sloppy paper had been published in a research journal where one 
finds many such examples. We have in fact discussed several of them in a 
series of communications [2-10]. In particular we want to draw the reader's 
attention to the extraordinary case of a predator-prey model that predicts a 
negative number of rabbits [5]. 

Finally, we mention that when a mild version of this article was submitted 
to the journal we were told " We urge you to contact the authors of this article 
before you submit a comment for publication. This approach is probably 
going to be much more productive for all concerned. I am going to reject the 
manuscript for now, but if after speaking directly with the authors you've 
come to an agreement that this manuscript should be published then you 
can resubmit." If we understand it clearly we are urged to ask permission 
from the authors to criticize their paper. Who would agree to it? 
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